[next] [prev] [prev-tail] [tail] [up]

\(\int \frac {\sqrt {\cos (c+d x)} (A+C \sec ^2(c+d x))}{(a+a \sec (c+d x))^{3/2}} \, dx\) [1163]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [B] (verified)
   Fricas [A] (verification not implemented)
   Sympy [F]
   Maxima [B] (verification not implemented)
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 37, antiderivative size = 172 \[ \int \frac {\sqrt {\cos (c+d x)} \left (A+C \sec ^2(c+d x)\right )}{(a+a \sec (c+d x))^{3/2}} \, dx=-\frac {(7 A-C) \text {arctanh}\left (\frac {\sqrt {a} \sqrt {\sec (c+d x)} \sin (c+d x)}{\sqrt {2} \sqrt {a+a \sec (c+d x)}}\right ) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}}{2 \sqrt {2} a^{3/2} d}-\frac {(A+C) \sin (c+d x)}{2 d \sqrt {\cos (c+d x)} (a+a \sec (c+d x))^{3/2}}+\frac {(5 A+C) \sin (c+d x)}{2 a d \sqrt {\cos (c+d x)} \sqrt {a+a \sec (c+d x)}} \]

[Out]

-1/2*(A+C)*sin(d*x+c)/d/(a+a*sec(d*x+c))^(3/2)/cos(d*x+c)^(1/2)-1/4*(7*A-C)*arctanh(1/2*sin(d*x+c)*a^(1/2)*sec
(d*x+c)^(1/2)*2^(1/2)/(a+a*sec(d*x+c))^(1/2))*cos(d*x+c)^(1/2)*sec(d*x+c)^(1/2)/a^(3/2)/d*2^(1/2)+1/2*(5*A+C)*
sin(d*x+c)/a/d/cos(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(1/2)

Rubi [A] (verified)

Time = 0.58 (sec) , antiderivative size = 172, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.135, Rules used = {4350, 4170, 4098, 3893, 212} \[ \int \frac {\sqrt {\cos (c+d x)} \left (A+C \sec ^2(c+d x)\right )}{(a+a \sec (c+d x))^{3/2}} \, dx=-\frac {(7 A-C) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)} \text {arctanh}\left (\frac {\sqrt {a} \sin (c+d x) \sqrt {\sec (c+d x)}}{\sqrt {2} \sqrt {a \sec (c+d x)+a}}\right )}{2 \sqrt {2} a^{3/2} d}+\frac {(5 A+C) \sin (c+d x)}{2 a d \sqrt {\cos (c+d x)} \sqrt {a \sec (c+d x)+a}}-\frac {(A+C) \sin (c+d x)}{2 d \sqrt {\cos (c+d x)} (a \sec (c+d x)+a)^{3/2}} \]

[In]

Int[(Sqrt[Cos[c + d*x]]*(A + C*Sec[c + d*x]^2))/(a + a*Sec[c + d*x])^(3/2),x]

[Out]

-1/2*((7*A - C)*ArcTanh[(Sqrt[a]*Sqrt[Sec[c + d*x]]*Sin[c + d*x])/(Sqrt[2]*Sqrt[a + a*Sec[c + d*x]])]*Sqrt[Cos
[c + d*x]]*Sqrt[Sec[c + d*x]])/(Sqrt[2]*a^(3/2)*d) - ((A + C)*Sin[c + d*x])/(2*d*Sqrt[Cos[c + d*x]]*(a + a*Sec
[c + d*x])^(3/2)) + ((5*A + C)*Sin[c + d*x])/(2*a*d*Sqrt[Cos[c + d*x]]*Sqrt[a + a*Sec[c + d*x]])

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
 /; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 3893

Int[Sqrt[csc[(e_.) + (f_.)*(x_)]*(d_.)]/Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)], x_Symbol] :> Dist[-2*b*(d/
(a*f)), Subst[Int[1/(2*b - d*x^2), x], x, b*(Cot[e + f*x]/(Sqrt[a + b*Csc[e + f*x]]*Sqrt[d*Csc[e + f*x]]))], x
] /; FreeQ[{a, b, d, e, f}, x] && EqQ[a^2 - b^2, 0]

Rule 4098

Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*
(B_.) + (A_)), x_Symbol] :> Simp[A*Cot[e + f*x]*(a + b*Csc[e + f*x])^m*((d*Csc[e + f*x])^n/(f*n)), x] - Dist[(
a*A*m - b*B*n)/(b*d*n), Int[(a + b*Csc[e + f*x])^m*(d*Csc[e + f*x])^(n + 1), x], x] /; FreeQ[{a, b, d, e, f, A
, B, m, n}, x] && NeQ[A*b - a*B, 0] && EqQ[a^2 - b^2, 0] && EqQ[m + n + 1, 0] &&  !LeQ[m, -1]

Rule 4170

Int[((A_.) + csc[(e_.) + (f_.)*(x_)]^2*(C_.))*(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b
_.) + (a_))^(m_), x_Symbol] :> Simp[(-a)*(A + C)*Cot[e + f*x]*(a + b*Csc[e + f*x])^m*((d*Csc[e + f*x])^n/(a*f*
(2*m + 1))), x] + Dist[1/(a*b*(2*m + 1)), Int[(a + b*Csc[e + f*x])^(m + 1)*(d*Csc[e + f*x])^n*Simp[b*C*n + A*b
*(2*m + n + 1) - (a*(A*(m + n + 1) - C*(m - n)))*Csc[e + f*x], x], x], x] /; FreeQ[{a, b, d, e, f, A, C, n}, x
] && EqQ[a^2 - b^2, 0] && LtQ[m, -2^(-1)]

Rule 4350

Int[(cos[(a_.) + (b_.)*(x_)]*(c_.))^(m_.)*(u_), x_Symbol] :> Dist[(c*Cos[a + b*x])^m*(c*Sec[a + b*x])^m, Int[A
ctivateTrig[u]/(c*Sec[a + b*x])^m, x], x] /; FreeQ[{a, b, c, m}, x] &&  !IntegerQ[m] && KnownSecantIntegrandQ[
u, x]

Rubi steps \begin{align*} \text {integral}& = \left (\sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \int \frac {A+C \sec ^2(c+d x)}{\sqrt {\sec (c+d x)} (a+a \sec (c+d x))^{3/2}} \, dx \\ & = -\frac {(A+C) \sin (c+d x)}{2 d \sqrt {\cos (c+d x)} (a+a \sec (c+d x))^{3/2}}-\frac {\left (\sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \int \frac {-\frac {1}{2} a (5 A+C)+a (A-C) \sec (c+d x)}{\sqrt {\sec (c+d x)} \sqrt {a+a \sec (c+d x)}} \, dx}{2 a^2} \\ & = -\frac {(A+C) \sin (c+d x)}{2 d \sqrt {\cos (c+d x)} (a+a \sec (c+d x))^{3/2}}+\frac {(5 A+C) \sin (c+d x)}{2 a d \sqrt {\cos (c+d x)} \sqrt {a+a \sec (c+d x)}}-\frac {\left ((7 A-C) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \int \frac {\sqrt {\sec (c+d x)}}{\sqrt {a+a \sec (c+d x)}} \, dx}{4 a} \\ & = -\frac {(A+C) \sin (c+d x)}{2 d \sqrt {\cos (c+d x)} (a+a \sec (c+d x))^{3/2}}+\frac {(5 A+C) \sin (c+d x)}{2 a d \sqrt {\cos (c+d x)} \sqrt {a+a \sec (c+d x)}}+\frac {\left ((7 A-C) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}\right ) \text {Subst}\left (\int \frac {1}{2 a-x^2} \, dx,x,-\frac {a \sqrt {\sec (c+d x)} \sin (c+d x)}{\sqrt {a+a \sec (c+d x)}}\right )}{2 a d} \\ & = -\frac {(7 A-C) \text {arctanh}\left (\frac {\sqrt {a} \sqrt {\sec (c+d x)} \sin (c+d x)}{\sqrt {2} \sqrt {a+a \sec (c+d x)}}\right ) \sqrt {\cos (c+d x)} \sqrt {\sec (c+d x)}}{2 \sqrt {2} a^{3/2} d}-\frac {(A+C) \sin (c+d x)}{2 d \sqrt {\cos (c+d x)} (a+a \sec (c+d x))^{3/2}}+\frac {(5 A+C) \sin (c+d x)}{2 a d \sqrt {\cos (c+d x)} \sqrt {a+a \sec (c+d x)}} \\ \end{align*}

Mathematica [A] (verified)

Time = 5.85 (sec) , antiderivative size = 339, normalized size of antiderivative = 1.97 \[ \int \frac {\sqrt {\cos (c+d x)} \left (A+C \sec ^2(c+d x)\right )}{(a+a \sec (c+d x))^{3/2}} \, dx=\frac {\sqrt {\sec (c+d x)} \left (7 \sqrt {2} A \arctan \left (\frac {\sqrt {2} \sqrt {\sec (c+d x)}}{\sqrt {1-\sec (c+d x)}}\right )-\sqrt {2} C \arctan \left (\frac {\sqrt {2} \sqrt {\sec (c+d x)}}{\sqrt {1-\sec (c+d x)}}\right )+\frac {8 A \sqrt {1-\sec (c+d x)}}{\sqrt {\sec (c+d x)}}+7 \sqrt {2} A \arctan \left (\frac {\sqrt {2} \sqrt {\sec (c+d x)}}{\sqrt {1-\sec (c+d x)}}\right ) \sec (c+d x)-\sqrt {2} C \arctan \left (\frac {\sqrt {2} \sqrt {\sec (c+d x)}}{\sqrt {1-\sec (c+d x)}}\right ) \sec (c+d x)+10 A \sqrt {-((-1+\sec (c+d x)) \sec (c+d x))}+2 C \sqrt {-((-1+\sec (c+d x)) \sec (c+d x))}+2 C \arcsin \left (\sqrt {1-\sec (c+d x)}\right ) (1+\sec (c+d x))+2 C \arcsin \left (\sqrt {\sec (c+d x)}\right ) (1+\sec (c+d x))\right ) \sin (c+d x)}{4 d \sqrt {-1+\cos (c+d x)} (a (1+\sec (c+d x)))^{3/2}} \]

[In]

Integrate[(Sqrt[Cos[c + d*x]]*(A + C*Sec[c + d*x]^2))/(a + a*Sec[c + d*x])^(3/2),x]

[Out]

(Sqrt[Sec[c + d*x]]*(7*Sqrt[2]*A*ArcTan[(Sqrt[2]*Sqrt[Sec[c + d*x]])/Sqrt[1 - Sec[c + d*x]]] - Sqrt[2]*C*ArcTa
n[(Sqrt[2]*Sqrt[Sec[c + d*x]])/Sqrt[1 - Sec[c + d*x]]] + (8*A*Sqrt[1 - Sec[c + d*x]])/Sqrt[Sec[c + d*x]] + 7*S
qrt[2]*A*ArcTan[(Sqrt[2]*Sqrt[Sec[c + d*x]])/Sqrt[1 - Sec[c + d*x]]]*Sec[c + d*x] - Sqrt[2]*C*ArcTan[(Sqrt[2]*
Sqrt[Sec[c + d*x]])/Sqrt[1 - Sec[c + d*x]]]*Sec[c + d*x] + 10*A*Sqrt[-((-1 + Sec[c + d*x])*Sec[c + d*x])] + 2*
C*Sqrt[-((-1 + Sec[c + d*x])*Sec[c + d*x])] + 2*C*ArcSin[Sqrt[1 - Sec[c + d*x]]]*(1 + Sec[c + d*x]) + 2*C*ArcS
in[Sqrt[Sec[c + d*x]]]*(1 + Sec[c + d*x]))*Sin[c + d*x])/(4*d*Sqrt[-1 + Cos[c + d*x]]*(a*(1 + Sec[c + d*x]))^(
3/2))

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(308\) vs. \(2(143)=286\).

Time = 0.76 (sec) , antiderivative size = 309, normalized size of antiderivative = 1.80

method result size
default \(\frac {\left (7 A \cos \left (d x +c \right ) \sqrt {2}\, \arctan \left (\frac {\sin \left (d x +c \right ) \sqrt {2}}{2 \left (1+\cos \left (d x +c \right )\right ) \sqrt {-\frac {1}{1+\cos \left (d x +c \right )}}}\right )+8 A \cos \left (d x +c \right ) \sin \left (d x +c \right ) \sqrt {-\frac {1}{1+\cos \left (d x +c \right )}}-C \cos \left (d x +c \right ) \sqrt {2}\, \arctan \left (\frac {\sin \left (d x +c \right ) \sqrt {2}}{2 \left (1+\cos \left (d x +c \right )\right ) \sqrt {-\frac {1}{1+\cos \left (d x +c \right )}}}\right )+7 A \sqrt {2}\, \arctan \left (\frac {\sin \left (d x +c \right ) \sqrt {2}}{2 \left (1+\cos \left (d x +c \right )\right ) \sqrt {-\frac {1}{1+\cos \left (d x +c \right )}}}\right )+10 A \sin \left (d x +c \right ) \sqrt {-\frac {1}{1+\cos \left (d x +c \right )}}-C \sqrt {2}\, \arctan \left (\frac {\sin \left (d x +c \right ) \sqrt {2}}{2 \left (1+\cos \left (d x +c \right )\right ) \sqrt {-\frac {1}{1+\cos \left (d x +c \right )}}}\right )+2 C \sin \left (d x +c \right ) \sqrt {-\frac {1}{1+\cos \left (d x +c \right )}}\right ) \sqrt {\cos \left (d x +c \right )}\, \sqrt {a \left (1+\sec \left (d x +c \right )\right )}}{4 a^{2} d \left (1+\cos \left (d x +c \right )\right )^{2} \sqrt {-\frac {1}{1+\cos \left (d x +c \right )}}}\) \(309\)

[In]

int((A+C*sec(d*x+c)^2)*cos(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(3/2),x,method=_RETURNVERBOSE)

[Out]

1/4/a^2/d*(7*A*cos(d*x+c)*2^(1/2)*arctan(1/2*sin(d*x+c)*2^(1/2)/(1+cos(d*x+c))/(-1/(1+cos(d*x+c)))^(1/2))+8*A*
cos(d*x+c)*sin(d*x+c)*(-1/(1+cos(d*x+c)))^(1/2)-C*cos(d*x+c)*2^(1/2)*arctan(1/2*sin(d*x+c)*2^(1/2)/(1+cos(d*x+
c))/(-1/(1+cos(d*x+c)))^(1/2))+7*A*2^(1/2)*arctan(1/2*sin(d*x+c)*2^(1/2)/(1+cos(d*x+c))/(-1/(1+cos(d*x+c)))^(1
/2))+10*A*sin(d*x+c)*(-1/(1+cos(d*x+c)))^(1/2)-C*2^(1/2)*arctan(1/2*sin(d*x+c)*2^(1/2)/(1+cos(d*x+c))/(-1/(1+c
os(d*x+c)))^(1/2))+2*C*sin(d*x+c)*(-1/(1+cos(d*x+c)))^(1/2))*cos(d*x+c)^(1/2)*(a*(1+sec(d*x+c)))^(1/2)/(1+cos(
d*x+c))^2/(-1/(1+cos(d*x+c)))^(1/2)

Fricas [A] (verification not implemented)

none

Time = 0.29 (sec) , antiderivative size = 406, normalized size of antiderivative = 2.36 \[ \int \frac {\sqrt {\cos (c+d x)} \left (A+C \sec ^2(c+d x)\right )}{(a+a \sec (c+d x))^{3/2}} \, dx=\left [-\frac {\sqrt {2} {\left ({\left (7 \, A - C\right )} \cos \left (d x + c\right )^{2} + 2 \, {\left (7 \, A - C\right )} \cos \left (d x + c\right ) + 7 \, A - C\right )} \sqrt {a} \log \left (-\frac {a \cos \left (d x + c\right )^{2} - 2 \, \sqrt {2} \sqrt {a} \sqrt {\frac {a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sqrt {\cos \left (d x + c\right )} \sin \left (d x + c\right ) - 2 \, a \cos \left (d x + c\right ) - 3 \, a}{\cos \left (d x + c\right )^{2} + 2 \, \cos \left (d x + c\right ) + 1}\right ) - 4 \, {\left (4 \, A \cos \left (d x + c\right ) + 5 \, A + C\right )} \sqrt {\frac {a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sqrt {\cos \left (d x + c\right )} \sin \left (d x + c\right )}{8 \, {\left (a^{2} d \cos \left (d x + c\right )^{2} + 2 \, a^{2} d \cos \left (d x + c\right ) + a^{2} d\right )}}, \frac {\sqrt {2} {\left ({\left (7 \, A - C\right )} \cos \left (d x + c\right )^{2} + 2 \, {\left (7 \, A - C\right )} \cos \left (d x + c\right ) + 7 \, A - C\right )} \sqrt {-a} \arctan \left (\frac {\sqrt {2} \sqrt {-a} \sqrt {\frac {a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sqrt {\cos \left (d x + c\right )}}{a \sin \left (d x + c\right )}\right ) + 2 \, {\left (4 \, A \cos \left (d x + c\right ) + 5 \, A + C\right )} \sqrt {\frac {a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sqrt {\cos \left (d x + c\right )} \sin \left (d x + c\right )}{4 \, {\left (a^{2} d \cos \left (d x + c\right )^{2} + 2 \, a^{2} d \cos \left (d x + c\right ) + a^{2} d\right )}}\right ] \]

[In]

integrate((A+C*sec(d*x+c)^2)*cos(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(3/2),x, algorithm="fricas")

[Out]

[-1/8*(sqrt(2)*((7*A - C)*cos(d*x + c)^2 + 2*(7*A - C)*cos(d*x + c) + 7*A - C)*sqrt(a)*log(-(a*cos(d*x + c)^2
- 2*sqrt(2)*sqrt(a)*sqrt((a*cos(d*x + c) + a)/cos(d*x + c))*sqrt(cos(d*x + c))*sin(d*x + c) - 2*a*cos(d*x + c)
 - 3*a)/(cos(d*x + c)^2 + 2*cos(d*x + c) + 1)) - 4*(4*A*cos(d*x + c) + 5*A + C)*sqrt((a*cos(d*x + c) + a)/cos(
d*x + c))*sqrt(cos(d*x + c))*sin(d*x + c))/(a^2*d*cos(d*x + c)^2 + 2*a^2*d*cos(d*x + c) + a^2*d), 1/4*(sqrt(2)
*((7*A - C)*cos(d*x + c)^2 + 2*(7*A - C)*cos(d*x + c) + 7*A - C)*sqrt(-a)*arctan(sqrt(2)*sqrt(-a)*sqrt((a*cos(
d*x + c) + a)/cos(d*x + c))*sqrt(cos(d*x + c))/(a*sin(d*x + c))) + 2*(4*A*cos(d*x + c) + 5*A + C)*sqrt((a*cos(
d*x + c) + a)/cos(d*x + c))*sqrt(cos(d*x + c))*sin(d*x + c))/(a^2*d*cos(d*x + c)^2 + 2*a^2*d*cos(d*x + c) + a^
2*d)]

Sympy [F]

\[ \int \frac {\sqrt {\cos (c+d x)} \left (A+C \sec ^2(c+d x)\right )}{(a+a \sec (c+d x))^{3/2}} \, dx=\int \frac {\left (A + C \sec ^{2}{\left (c + d x \right )}\right ) \sqrt {\cos {\left (c + d x \right )}}}{\left (a \left (\sec {\left (c + d x \right )} + 1\right )\right )^{\frac {3}{2}}}\, dx \]

[In]

integrate((A+C*sec(d*x+c)**2)*cos(d*x+c)**(1/2)/(a+a*sec(d*x+c))**(3/2),x)

[Out]

Integral((A + C*sec(c + d*x)**2)*sqrt(cos(c + d*x))/(a*(sec(c + d*x) + 1))**(3/2), x)

Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 8312 vs. \(2 (143) = 286\).

Time = 0.57 (sec) , antiderivative size = 8312, normalized size of antiderivative = 48.33 \[ \int \frac {\sqrt {\cos (c+d x)} \left (A+C \sec ^2(c+d x)\right )}{(a+a \sec (c+d x))^{3/2}} \, dx=\text {Too large to display} \]

[In]

integrate((A+C*sec(d*x+c)^2)*cos(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(3/2),x, algorithm="maxima")

[Out]

-1/4*((4*(7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*
d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x +
3/2*c)^4 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2
*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^4 + 4*(7*log(cos(
1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*d*x + 1/2*c)^2 + sin
(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c)^4 + 70*(log(c
os(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + si
n(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^2*sin(1/2*d*x + 1/2*c)^2 + 7*(log(cos
(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(
1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^4 - 8*sin(1/2*d*x + 1/2*c)^5 + 28*(7*(l
og(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2
+ sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*
d*x + 1/2*c))*cos(3/2*d*x + 3/2*c)^3 + 4*(21*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*
d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2
*d*x + 1/2*c) - 24*sin(1/2*d*x + 1/2*c)^2 - 20)*sin(3/2*d*x + 3/2*c)^3 - 8*(10*cos(1/2*d*x + 1/2*c)^2 + 3)*sin
(1/2*d*x + 1/2*c)^3 + ((7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) -
7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*
cos(3/2*d*x + 3/2*c)^2 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1)
 - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^2 +
 (7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*d*x + 1/
2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c)^2
 + 7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/
2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2 - 8*sin(1/2*d*x + 1/2*c)
^3 + 6*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x
 + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2
*c)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c) + 2*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2
*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1)
)*sin(1/2*d*x + 1/2*c) - 8*sin(1/2*d*x + 1/2*c)^2 - 8)*sin(3/2*d*x + 3/2*c) - 8*(9*cos(1/2*d*x + 1/2*c)^2 + 2)
*sin(1/2*d*x + 1/2*c))*cos(5/2*d*x + 5/2*c)^2 + (427*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*
sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))
*cos(1/2*d*x + 1/2*c)^2 + 35*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1
) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2
- 40*sin(1/2*d*x + 1/2*c)^3 - 8*(61*cos(1/2*d*x + 1/2*c)^2 + 9)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c)^2 +
 ((7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*d*x + 1
/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c)^
2 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x +
1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^2 + (7*log(cos(1/2*d*x +
 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x
+ 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c)^2 + 7*(log(cos(1/2*d*x
 + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x
+ 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2 - 8*sin(1/2*d*x + 1/2*c)^3 + 6*(7*(log(cos(1/
2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2
*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2
*c))*cos(3/2*d*x + 3/2*c) + 2*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c)
 + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c
) - 8*sin(1/2*d*x + 1/2*c)^2 - 8)*sin(3/2*d*x + 3/2*c) - 8*(9*cos(1/2*d*x + 1/2*c)^2 + 2)*sin(1/2*d*x + 1/2*c)
)*sin(5/2*d*x + 5/2*c)^2 + (8*(7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c)
+ 1) - 7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1
/2*c))*cos(3/2*d*x + 3/2*c)^2 + 259*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2
*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/
2*c)^2 + 91*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d
*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2 - 104*sin(1/2*d*x
 + 1/2*c)^3 + 28*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(c
os(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2
*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c) - 8*(37*cos(1/2*d*x + 1/2*c)^2 + 21)*sin(1/2*d*x + 1/
2*c))*sin(3/2*d*x + 3/2*c)^2 + 2*(2*(7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1
/2*c) + 1) - 7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d
*x + 1/2*c))*cos(3/2*d*x + 3/2*c)^3 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x
+ 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x
 + 1/2*c)^3 + 7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1
/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1
/2*c)^2 - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^3 + 13*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/
2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1
/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c)^2 + (2*(7
*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*d*x + 1/2*c
)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c) + 7*
(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^
2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/
2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c)^2 + 2*(84*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/
2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1
/2*d*x + 1/2*c)^2 + 7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log
(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2 - 8*sin
(1/2*d*x + 1/2*c)^3 - 16*(6*cos(1/2*d*x + 1/2*c)^2 + 1)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c) + 2*(7*(log
(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 +
sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c) - 8*cos(1/2*d*
x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2 + 2*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x
+ 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x
 + 1/2*c) - 8*sin(1/2*d*x + 1/2*c)^2 - 8)*cos(3/2*d*x + 3/2*c) - 8*cos(1/2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c)
- 8*(9*cos(1/2*d*x + 1/2*c)^3 + 2*cos(1/2*d*x + 1/2*c))*sin(1/2*d*x + 1/2*c))*cos(5/2*d*x + 5/2*c) + 2*(147*(l
og(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2
+ sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^3 + 35*(log(cos(1/2*d*x + 1/2*c)^
2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2
 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2 - 40*cos(1/2*d*x + 1/2*c)*sin(1/2*
d*x + 1/2*c)^3 - 56*(3*cos(1/2*d*x + 1/2*c)^3 + cos(1/2*d*x + 1/2*c))*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*
c) + 2*(2*(7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2
*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*sin(3/2*d*x +
 3/2*c)^3 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/
2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^2*sin(1/2*d*x +
1/2*c) + 7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*
x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^3 - 8*sin(1/2*d*x +
1/2*c)^4 + (7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2
*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c) - 8*sin(1/2*d*x +
 1/2*c)^2 - 4)*cos(3/2*d*x + 3/2*c)^2 + (35*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d
*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*
d*x + 1/2*c) - 40*sin(1/2*d*x + 1/2*c)^2 - 36)*sin(3/2*d*x + 3/2*c)^2 - 4*(18*cos(1/2*d*x + 1/2*c)^2 + 5)*sin(
1/2*d*x + 1/2*c)^2 + 6*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) -
 log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)*sin(1
/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2 - 4*cos(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c)
 - 36*cos(1/2*d*x + 1/2*c)^2 + 2*((7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2
*c) + 1) - 7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x
 + 1/2*c))*cos(3/2*d*x + 3/2*c)^2 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x +
1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x +
 1/2*c)^2 + 14*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/
2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2 - 16*sin(1/2*d
*x + 1/2*c)^3 + 6*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(
cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/
2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c) - 4*(18*cos(1/2*d*x + 1/2*c)^2 + 7)*sin(1/2*d*x + 1/
2*c))*sin(3/2*d*x + 3/2*c))*sin(5/2*d*x + 5/2*c) + 2*(133*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2
 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c)
+ 1))*cos(1/2*d*x + 1/2*c)^2*sin(1/2*d*x + 1/2*c) + 21*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 +
2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1
))*sin(1/2*d*x + 1/2*c)^3 - 24*sin(1/2*d*x + 1/2*c)^4 + 2*(21*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*
c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2
*c) + 1))*sin(1/2*d*x + 1/2*c) - 24*sin(1/2*d*x + 1/2*c)^2 - 20)*cos(3/2*d*x + 3/2*c)^2 - 8*(19*cos(1/2*d*x +
1/2*c)^2 + 7)*sin(1/2*d*x + 1/2*c)^2 + 16*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*
d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2
*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2 - 5*cos(1/2*d*x + 1/2*c))*c
os(3/2*d*x + 3/2*c) - 80*cos(1/2*d*x + 1/2*c)^2)*sin(3/2*d*x + 3/2*c) - 8*(9*cos(1/2*d*x + 1/2*c)^4 + 11*cos(1
/2*d*x + 1/2*c)^2)*sin(1/2*d*x + 1/2*c))*A*sqrt(a)/(4*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^4 + 28*sqrt(2)*a^2*cos(
3/2*d*x + 3/2*c)^3*cos(1/2*d*x + 1/2*c) + 9*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^4 + 4*sqrt(2)*a^2*sin(3/2*d*x + 3
/2*c)^4 + 12*sqrt(2)*a^2*sin(3/2*d*x + 3/2*c)^3*sin(1/2*d*x + 1/2*c) + 10*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2*s
in(1/2*d*x + 1/2*c)^2 + sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^4 + (sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^2 + 6*sqrt(2)*a
^2*cos(3/2*d*x + 3/2*c)*cos(1/2*d*x + 1/2*c) + 9*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2 + sqrt(2)*a^2*sin(3/2*d*x
+ 3/2*c)^2 + 2*sqrt(2)*a^2*sin(3/2*d*x + 3/2*c)*sin(1/2*d*x + 1/2*c) + sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^2)*cos
(5/2*d*x + 5/2*c)^2 + (61*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2 + 5*sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^2)*cos(3/2*d
*x + 3/2*c)^2 + (sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^2 + 6*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)*cos(1/2*d*x + 1/2*c)
+ 9*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2 + sqrt(2)*a^2*sin(3/2*d*x + 3/2*c)^2 + 2*sqrt(2)*a^2*sin(3/2*d*x + 3/2*
c)*sin(1/2*d*x + 1/2*c) + sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^2)*sin(5/2*d*x + 5/2*c)^2 + (8*sqrt(2)*a^2*cos(3/2*
d*x + 3/2*c)^2 + 28*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)*cos(1/2*d*x + 1/2*c) + 37*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c
)^2 + 13*sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^2)*sin(3/2*d*x + 3/2*c)^2 + 2*(2*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^3
+ 13*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^2*cos(1/2*d*x + 1/2*c) + 9*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^3 + sqrt(2)*
a^2*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2 + (2*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c) + sqrt(2)*a^2*cos(1/2*d*
x + 1/2*c))*sin(3/2*d*x + 3/2*c)^2 + 2*(12*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2 + sqrt(2)*a^2*sin(1/2*d*x + 1/2*
c)^2)*cos(3/2*d*x + 3/2*c) + 2*(2*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)*sin(1/2*d*x + 1/2*c) + sqrt(2)*a^2*cos(1/2*
d*x + 1/2*c)*sin(1/2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c))*cos(5/2*d*x + 5/2*c) + 2*(21*sqrt(2)*a^2*cos(1/2*d*x
+ 1/2*c)^3 + 5*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2)*cos(3/2*d*x + 3/2*c) + 2*(2*sqrt(2)*a^
2*sin(3/2*d*x + 3/2*c)^3 + sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^2*sin(1/2*d*x + 1/2*c) + 6*sqrt(2)*a^2*cos(3/2*d*x
 + 3/2*c)*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c) + 9*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2*sin(1/2*d*x + 1/2*c
) + 5*sqrt(2)*a^2*sin(3/2*d*x + 3/2*c)^2*sin(1/2*d*x + 1/2*c) + sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^3 + 2*(sqrt(2
)*a^2*cos(3/2*d*x + 3/2*c)^2 + 6*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)*cos(1/2*d*x + 1/2*c) + 9*sqrt(2)*a^2*cos(1/2
*d*x + 1/2*c)^2 + 2*sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^2)*sin(3/2*d*x + 3/2*c))*sin(5/2*d*x + 5/2*c) + 2*(6*sqrt
(2)*a^2*cos(3/2*d*x + 3/2*c)^2*sin(1/2*d*x + 1/2*c) + 16*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)*cos(1/2*d*x + 1/2*c)
*sin(1/2*d*x + 1/2*c) + 19*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2*sin(1/2*d*x + 1/2*c) + 3*sqrt(2)*a^2*sin(1/2*d*x
 + 1/2*c)^3)*sin(3/2*d*x + 3/2*c)) - (4*(sin(3/2*d*x + 3/2*c) - sin(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*
d*x + 3/2*c))))*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 8*(sin(3/2*d*x + 3/2*c) - sin(1
/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))))*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3
/2*c))) + (2*(2*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 1)*cos(4/3*arctan2(sin(3/2*d*x
+ 3/2*c), cos(3/2*d*x + 3/2*c))) + cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 4*cos(2/3*
arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3
/2*c)))^2 + 4*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))*sin(2/3*arctan2(sin(3/2*d*x + 3/2*c
), cos(3/2*d*x + 3/2*c))) + 4*sin(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 4*cos(2/3*arcta
n2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 1)*log(cos(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3
/2*c)))^2 + sin(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 2*sin(1/3*arctan2(sin(3/2*d*x + 3
/2*c), cos(3/2*d*x + 3/2*c))) + 1) - (2*(2*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 1)*c
os(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*
x + 3/2*c)))^2 + 4*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + sin(4/3*arctan2(sin(3/2*d*
x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 4*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))*sin(2/3*
arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 4*sin(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3
/2*c)))^2 + 4*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 1)*log(cos(1/3*arctan2(sin(3/2*d*
x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + sin(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 - 2*sin(
1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 1) - 4*(cos(3/2*d*x + 3/2*c) - cos(1/3*arctan2(sin(
3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))))*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) - 8*(co
s(3/2*d*x + 3/2*c) - cos(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))))*sin(2/3*arctan2(sin(3/2*d*x
 + 3/2*c), cos(3/2*d*x + 3/2*c))) + 4*sin(3/2*d*x + 3/2*c) - 4*sin(1/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d
*x + 3/2*c))))*C/((sqrt(2)*a*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + 4*sqrt(2)*a*cos(
2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))^2 + sqrt(2)*a*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), c
os(3/2*d*x + 3/2*c)))^2 + 4*sqrt(2)*a*sin(4/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c)))*sin(2/3*arc
tan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 4*sqrt(2)*a*sin(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*
d*x + 3/2*c)))^2 + 4*sqrt(2)*a*cos(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + 2*(2*sqrt(2)*a*c
os(2/3*arctan2(sin(3/2*d*x + 3/2*c), cos(3/2*d*x + 3/2*c))) + sqrt(2)*a)*cos(4/3*arctan2(sin(3/2*d*x + 3/2*c),
 cos(3/2*d*x + 3/2*c))) + sqrt(2)*a)*sqrt(a)))/d

Giac [F]

\[ \int \frac {\sqrt {\cos (c+d x)} \left (A+C \sec ^2(c+d x)\right )}{(a+a \sec (c+d x))^{3/2}} \, dx=\int { \frac {{\left (C \sec \left (d x + c\right )^{2} + A\right )} \sqrt {\cos \left (d x + c\right )}}{{\left (a \sec \left (d x + c\right ) + a\right )}^{\frac {3}{2}}} \,d x } \]

[In]

integrate((A+C*sec(d*x+c)^2)*cos(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(3/2),x, algorithm="giac")

[Out]

integrate((C*sec(d*x + c)^2 + A)*sqrt(cos(d*x + c))/(a*sec(d*x + c) + a)^(3/2), x)

Mupad [F(-1)]

Timed out. \[ \int \frac {\sqrt {\cos (c+d x)} \left (A+C \sec ^2(c+d x)\right )}{(a+a \sec (c+d x))^{3/2}} \, dx=\int \frac {\sqrt {\cos \left (c+d\,x\right )}\,\left (A+\frac {C}{{\cos \left (c+d\,x\right )}^2}\right )}{{\left (a+\frac {a}{\cos \left (c+d\,x\right )}\right )}^{3/2}} \,d x \]

[In]

int((cos(c + d*x)^(1/2)*(A + C/cos(c + d*x)^2))/(a + a/cos(c + d*x))^(3/2),x)

[Out]

int((cos(c + d*x)^(1/2)*(A + C/cos(c + d*x)^2))/(a + a/cos(c + d*x))^(3/2), x)

[next] [prev] [prev-tail] [front] [up]